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1. Introduction 


In this paper we prove boundary regularity estimates for finite difference schemes for 
elliptic systems of partial differential equations. These estimates express the smoothness, 
or regularity, of the solution of the difference scheme up to the boundary in terms of the 
smoothness of the data in the interior and on the boundary of the domain on which the 
equations hold. For those difference schemes which require more boundary conditions than 
the corresponding differential equations, the estimates show the effect of the additional 
boundary data on the smoothness of the solution. 

Our results are analogous to the Sobolev norm estimates of Agmon, Douglis and 
Nirenberg [l] for solutions of elliptic systems of differential equations. In previous work 
(Bube and Strikwerda [3]) interior regularity estimates for systems of elliptic difference 
schemes were derived and a theory of pseudo— difference operators was developed (see also 
Frank [4]). It is assumed that the reader is familiar with this paper. That theory is used 
here to prove the boundary regularity estimates. 

The approach used in this paper is similar to that used by Gustafsson, Kreiss, and 
Sundstrom [6] on the initial boundary value problem for difference schemes for hyperbolic 
equations. The essential ideas are first to transform the system to a one-step scheme in 
the direction normal to the boundary, and then construct a Hermitian pseudo-difference 
operator, the Gustafsson-Kreiss-Sundstrom symmetrizer, for this one-step scheme which 
enables one to obtain the appropriate estimates for the solution. Michelson [9] has extended 
these ideas to multidimensional initial-boundary value problems. Because the principal 
symbol of an elliptic system of partial differential equations can involve different orders of 
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differentiation on different variables, the construction used here is more general than that 
of Gustafsson, Kreiss, and Sundstrom. 

We consider only boundary-fitted grids, i.e., those in which the boundary is a co- 
ordinate surface. This excludes the grid systems used by Bramble and Hubbard [2] and 
others in which the boundary curve is not parallel to a coordinate line. Boundary-fitted 
grids are in common use in computational fluid dynamics (e.g. Thompson et al. [12]). 
A scheme for the Stokes equations, an elliptic system, which is second-order accurate on 
boundary-fitted grids has been presented by Strikwerda [10], 

The significance of the regularity estimates proven here is that they show the effect of 
the discretizations of the boundary conditions on the smoothness of the solution near the 
boundary. Regular elliptic finite difference schemes satisfy interior regularity estimates 
and thus have smooth solutions away from the boundary. In this paper the conditions 
under which the solutions have optimal smoothness at the boundary are given. 

In a subsequent work we intend to extend the results given here to estimates of the 
accuracy of the finite difference solution as an approximation to the solution of an elliptic 
system of partial differential equations. In this paper we give only the estimate in Theorem 
3.2 which gives the accuracy of the divided differences of the solution in terms of the 
accuracy of the solution itself. 

The outline of the paper is as follows. We begin by considering elliptic systems of 
difference schemes defined on a half— space with a rectangular grid. Section 2 is concerned 
with the basic definitions and assumptions. The main result of this paper, Theorem 
3.1, is stated in Section 3 after the Complementing Conditions are defined. Section 4 is 
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concerned with the reduction of general schemes to a one-step scheme, and in Section 5 
the Gustafsson-Kreiss-Sundstrom symmetrizer is constructed. The boundary regularity 
estimates are obtained in Section 6 which proves the main theorem, Theorem 3.1. Section 
7 is a discussion of the more general situation of difference schemes on domains with 
curvilinear grids, and Section 8 is a summary of the paper which is intended to be a user’s 
guide, enabling one to apply the results of the paper to particular problems. 

2. Definitions and Assumptions 

We begin by considering the regular elliptic system of difference equations 

n 

^ ^ LijUj (xj,, y n) — Fi^Xf/, i l,...,n, n £ Z , (2-1) 

y=i 

on the half-space R^. = {(x,y) : x > 0,y G R d_1 ,d > 2}, where the grid points are 
x„ = uh, = fih, v > 0. The boundary conditions are given by 

n 

^ BkjUj (-^o> Vfi) = *ftfc(i//t)> fc = l,...,9j € Z , (2*2) 

i=i 

where the Bkj are difference operators which involve only forward differences and trans- 
lations. The number of boundary conditions which must be specified will be determined 
later, (see Assumption 2.1). Without loss of generality (see Section 7), we assume that 
F(x,y) and <j>k{y) are -periodic in each yj. The system (2.1) is assumed to be regular 
elliptic of order (<r, r), (Bube and Strikwerda [3]). We also assume that there is a p G Z 9 
such that 

e s " +r '- < 2 - 3 ) 
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where the b^j are the symbols of the B^j when considered as difference operators on all of 

R d . 

Following Bube and Strikwerda [3], we define the tangential grid 
G •“ { J//x : 0 ^ l^j < 27V , j = 1 , d — 1}, 
where := h/i G [0, 2tt] 1 and h := 2tt/( 27V + 1). We also define 

(«(x,-),v(x,.)) G : = {^) d ~ x E u ( x >y)* v ( x >y)> 

y€G 

and 

y€G 

where r := {£ £ Z d ~ l : |£y| < N for j = 1 l}. By the Fourier inversion formula, 

u (x,y) = E e * y *u(x,$). 

ter 

The norms used in the regularity estimate are as follows. Let r € Z d_1 and s £ R. Define 

d-i 4 , 

A(/i,w) 2 := 1 + E — sin 2 (-w y ), u = h£ £ hT, 
i=i 

and 

IMS - £ £ £ £ |A(k,fc{)*- m c + 4i(*,.«p, 

i /=0 m =0 

ll«HJ+r :=ElKII^- 

y=i 
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At the boundary we define 


u - 2 - 

U J s ’ 


2 ^ 2 ^ 

£er m=o 


l 6™ + Uj(xo,t)\' 


l U ls+r l tt jla+ry‘ 

i - 1 

Note that |uy| s is the sum of norms of uy and its forward divided differences in x of order 
at most s evaluated at x 0 , and considered as discrete functions of y. 

We will assume that the operators L,y and Bkj contain only differences of order oy+ry 
and pk + Tj, respectively, with coefficients which are independent of h. That is, /,y, the 
symbol of L,y, is in S £T<+r '» and 

lij(h,x,y,$) =l' ij (x,y,s)h~( a<+T i ) , (2.4) 

f = teZ d , |e| <N, h = 2n/{2N+l). 

for some function independent of h. Similarly for bkj{h, y, f). The inclusion of lower 
order terms does not affect the form of the final regularity estimates, as will be discussed 
in Section 7. For simplicity of exposition we will assume that the symbols lij(h,x,y,$) are 
independent of x. 

We will also assume that Tj > 1 for each j and that max, a,- = 0. Note that if, for 
some j, Tj = 0 then «y may be expressed as a linear combination of differences of the other 
components of u, therefore uy may be eliminated from the system without affecting the 
ellipticity. 
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To apply the theory of pseudo-difference operators to the boundary value problem 
(2.1-2.2) we introduce the reduced system of difference equations as follows. Each of the 
difference operators L {j is a difference operator in both x and y. The reduced operator 
Li](h,y,u)) is the symbol of with respect to y only, and thus remains as a difference 
operator in x. Here u = h£, where f G Z d-1 is the dual variable to y. Boundary operators 
B kj{h,y,u) are defined in a similar way. The reduced system is 

n 

^2Lij{h,y,u})uj{x u , u) = Fi{x u ,u)), i = 1 , ( 2 . 5 ) 
;=i 

with boundary conditions 

n 

Yl^kj{h,y,u)uj(x 0 ,u) = <j>k{u), k = 1 (2.6) 

3 = 1 

The reduced system (2.5) will be transformed to a one-step difference scheme in x 
using a new vector of variables iy(x^,u;) so that the reduced system is equivalent to one 
of the form 


^(xj/+i,w) = M[y,oj)W(x l/ ,u) + hT(x l/ ,oj), v > 0, (2.7) 

with boundary conditions 

0(y,w)W(x o ,u>) = $(w). (2.8) 

The regularity estimates will be derived for the one-step scheme (2.7) using the func- 
tion W(x v , u). The final form of the estimate will follow by transferring back to the original 
function u(x„,u). 
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The estimates will be derived through the use of a matrix symbol H(y,u), the sym- 


metrizer, which will be constructed to satisfy the two matrix inequalities 


M*HM - H >c 0 h A 0 

(2.9) 

H + c^B > c 2 

(2.10) 


for some positive constants co,ci,C 2 and Ao = Ao(h,u>) := \j4h~ 2 s in 2 (|wi). The 
construction of the matrix H is analogous to that used by Gustafsson, Kreiss, and Sund- 
strom [6] for hyperbolic difference schemes. 

Before transforming the reduced equation to the one-step scheme we require some ad- 
ditional definitions and assumptions. Related to the reduced equation (2.5) is the resolvent 
equation 

n 

^2Uj{h,y,oj,z)vj = 0, i = l,...,n (2.11) 

i=i 

where /~,y(h,y,u;, z) = /,y(A, y, f), the symbol of L t; , and f = (} log(z),u>), Vj € C,j = 
1 ,...,n. Here log(z) is any logarithm of z €. C. 

Definition 2.1 

The values of z = z(y,u>) for which the resolvent equation (2.11) has nontrivial solu- 
tions are called eigenvalues of the resolvent equation. 

Let 

n 

R(y,u>,z) := det{lij(h,y,u,z)}h 2p , 2 p := ^(r,- + <r,). 

*=i 

Note that R(y,u,z) is independent of h by equation (2.4). The eigenvalues z(y,u) are the 
roots of R(y,u),z) = 0. We will need the following lemma. 
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Lemma 2.1 


The eigenvalues z(y,u>) satisfy 

l z (y> u, )| 7 ^ 1 if w^O. 

Moreover, if u> = 0 and \z\ — 1 then z — 1. 

Proof: 

Suppose z(y , cj ) = exp^'ft) with Cl real and |frj < n. Then 

0 = R(y,u,exp (*Vi)) = det{/,y(/i,y,w,exp (tfi))} 

implies 

det{/ tJ (/i,y,f)} = 0, 

where f = Since the matrix of the /, 7 is a regular elliptic symbol with only highest 

order terms, the determinant can not vanish unless C = 0. This proves the lemma. 

To transform the reduced equation (2.5) to the one-step scheme (2.7) we must place 
restrictions on the size of the stencil of the difference equations. 

Definition 2.2 

If Q is a difference operator in x written in the form 

b 

Qf{x„) = Y,Q^x)T»f( Xl ,), 

H~CL 

where neither Q a nor Q are identically zero, then the extent of Q is the ordered pair 
(a, 6 ). If Qi is a difference operator with extent ( ai,&! ) we say the extent of Qi is less 
than the extent of Q if a < aj < <6. 

The condition we place on the extent of the reduced system is: 
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Resolvent Condition 


If the operators Lij in the reduced equation (2.5) have extent (a,y,6 t y), we assume 
there are a - , 0 £ + ,/? - ,/? + in Z n such that 


jd j ^ &{j ^ ^ 0 ! t - *f* /5y , 1 5 ...,n 


( 2 . 12 ) 


and such that the number of roots z(y,oj), counting multiplicity, of the equation 


R(y,u,z) = 0 


is precisely 

J2( a i +Pk)-( a k +Pk) 

k = l 

for any value of w. If Lij is identically zero for some values of (t,j) then we place no 
restriction of the form (2.12) for that value of 

We now show that without loss of generality we may take a~ = (3~ = 0 in the 
resolvent condition. By subtracting some positive integer from all the and adding it to 
all the fdj we can have /3~ > 0, with min/?“ = 0, without altering the resolvent condition. 

By operating on the i-th equation with T h Q ‘ , where Th is the translation operator in 
the x direction, we obtain an equivalent system with a~ = 0. (The new value of o: + is 
— a - .) Then by defining new dependent variables 

u'j{x v ,uj) = Uj(x u + p~h,u>), v > 0 

a new system is obtained with /3~ = 0. (The new value of (3 + is /? + — f3~.) Note that 
the (old) variables uj{xo,u), ...,uy(i^-_ 1 ,w) do not appear in the difference equations and 
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thus are superfluous, appearing only in the boundary conditions. We will assume that these 
boundary variables can be expressed as linear combinations of the other (nonsuperfluous) 
variables and thus be eliminated. If these variables can not be so eliminated then the 
system does not have sufficient boundary conditions to determine the solution. 

We now make an assumption on the number of boundary conditions. 

Assumption 2.1 

The resolvent condition holds with a~ = /3~ = 0 and the number of boundary condi- 
tions q is equal to the number of roots in z of 


which satisfy 


R(y,u,z) = 0 


0 < |z(t/,u>)| < 1 for |w| 7^ 0. 


(2.13) 


Recall that \z{y,u)\ 7^ 1 for u ^ 0 by Lemma 2.1. We also know that 

Oi + Tj < af + (2.14) 

since a consistent difference operator approximating a differential operator of order s must 
involve at least s + 1 points. In the case that the number of boundary conditions q is larger 
than p as defined in Definition 2.1, we need an additional assumption. We assume that 
the boundary conditions are ordered so that the p k are in increasing order and we then 
define two important quantities. Let 

p := max (p k + 1,0). (2.15) 
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and 


f minfc> g (/)fc) + 1, if q > p; 
1 oo, if q = P- 


(2.16) 


The number p* is only used to limit the order of divided differences in the regularity 
estimates, when p* is infinite there is no restriction on the order of the differences, e.g. 
Theorem 3.1. 


Assumption 2.2 

If q > p, then we assume that 


Pk> P, for P < k < q. 


That is, the last q — p boundary conditions have weights pk which are nonnegative 
and larger than the weights of the first p boundary conditions. 

This assumption is needed to obtain the regularity estimates in the appropriate norms, 
and can be motivated as follows. If the difference scheme (2.1) is an approximation to a 
system of differential equations then p of the boundary conditions (2.2) would correspond 
to the boundary conditions of the differential equations. The remaining q — p boundary 
conditions which are required by the difference equations should be intrinsically distinct 
from the first p boundary conditions, and this distinction is maintained by having the 
weights pk of the last q — p boundary conditions sufficiently large and indeed larger than 
the weights of the first p boundary conditions. Note that p is zero for the classical Dirichlet 
and Neumann boundary value problems for elliptic equations of order 2p. 
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Definition 2.3 


Corresponding to the system of difference equations (2.1) is the associated system of 
differential equations 


22£ij(v,d x ,d v )wj(x,y) = Fi(x,y), t = l,...,n, (2.17) 

3 = 1 

A A 

where Lij is the differential operator whose symbol Uj(y, £) is the limit of lij(y, h£) as h 
tends to zero. The associated boundary conditions are 

n 

^B fc y(y,a 3 .| I=o ,5 y )u)_ ? (0,y) = ^ fc (y) k = l,...,p, (2.18) 

i=i 

These are obtained by the same limiting procedure from the first p boundary conditions 
(2.2), (see Assumption 2.2). 

The above limits exist by equation (2.4) since the Lij and bkj are difference operators 
in S <T,+T > and S Pt+7 > respectively; see also equation (2.8) of Bube and Strikwerda [3]. 

As in Agmon et al. [1, p.39], we require the associated differential equation to satisfy 
the following condition. 

Assumption 2.3 Supplementary condition on L. 

L{y, £) is of even degree 2 p with respect to £. For every pair of vectors £ and £' in R d 
the polynomial L(y,£ + t£') in the complex variable r has exactly p roots with positive 
imaginary part. 
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3. Complementing Condition 


In order for the regularity estimates to hold for the system (2.1) certain conditions 
must be satisfied by the boundary conditions (2.2). Before stating these conditions, called 
Complementing Conditions, it will be helpful to introduce some notation. We will write 
the reduced equation (2.5) as 


L(h, y,u))u(xi,,cj) = F(x u ,u), u > 0, 


(3.1) 


and the boundary conditions (2.6) as 


jBi(/i,y,w)ti(xo,w) = <£i(w), 
B 2 (h,y,u)u(x 0 ,u) = (foM, 


(3.2) 


where B i is composed of the first p boundary operators and Bi the last q — p boundary 
operators (see assumption 2.2). 

The reduced equation (3.1) can be replaced by an equivalent equation independent 
of h by the following scaling procedure. Multiply the i-th equation by h a ' and replace 
the variables uy by u’-h T i . This gives a system of difference equations in the u' which 
is equivalent to the original system. The boundary conditions can be scaled in a similar 
manner. 

The Complementing Conditions will be stated in terms of eigensolutions, of which 
there are three types. Because of the above scaling procedure, we need only consider h 
equal to 1. 
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Definition 3.1 


An eigensolution of type I is a nontrivial solution to the difference equation 
L(l,y,ut)u(x l/ ,uj) = 0 u>0, for some w^O, 


satisfying 

Bi{l,y,u)u(x 0 ,oj) = 0 

a) 

B 2 (l,y,u)u(x 0 ,u>) = 0, 

b) u(x v ,u>) -> 0 as v — » <x). 


Definition 3.2 

An eigensolution of type II is a nontrivial solution to the associated differential 

tion 

L{y,d x ,9)w(x,0) = 0 for 0 G R^ -1 , |0| = 1, 

satisfying 

a) B 1 (y,d x ,0)w(x,0)f x=o = 0, 

b) w(x, 0) —> 0 as x —> oo. 


In the case q > p we formulate: 

Definition 3.3 

An eigensolution of type III is a nontrivial solution to the difference equation 

y> |u=o = Oj 


equa- 
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satisfying 


а) 5 2 (l,y,xo)u(xo,w)| w=0 = 0, 

б) u(x v ,0) — > 0 as v — > oo. 


We now state the Complementing Condition and the main theorem of this paper. 
Complementing Condition 

The system (2.1) with boundary conditions (2.2) satisfies the Complementing Condi- 
tion if there are no eigensolutions of type I, II, or III. 

Theorem 3.1 

If u(x l/ ,y tl ) is a solution to the system (2.1) with the boundary conditions (2.2) and 
Assumptions 2.1, 2.2, and 2.3 are satisfied, then the following regularity estimate holds 
for each s with p < s < p* and h sufficiently small if, and only if, the Complementing 
Condition holds. 


u 


r-f s 


+ M r+s-i — CsG^ll 3-0-i + \h P t+ 2 < l> 2 \l-t + Ill’ll s—<r + IMIo)- 


(3.3) 


where t = p + f |_2(s — p) J. 

The boundary data (f>\ and <^ 2 are defined by equation (3.2) according to Assumption 2.2. 
The subscript on the norm of <j>\ uses only the first p components of the multi-index p, 
and h p ~ aJr *4>2 is understood as multiplication of the (k — p + l)-st component of <£ 2 by 
k pk ~ s+ 2 for k > p. 

Note that the nonexistence of eigensolutions of type II is equivalent to a regularity 
estimate analogous to (3.3) (with $ 2 = 0 ) holding for the associated differential equation 
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(2.14). Theorem 3.1 is proved in Section 6 as Theorems 6.1 and 6.2. 

An immediate corollary of Theorem 3.1 is Theorem 3.2. 

Theorem 3.2 

Suppose v is a solution to the finite difference equations (2.1) and u is a solution to 
the corresponding partial differential equation (2.17) such that 

1) Lv = F, Bv = <j>, 

2) Lu = F+0{h r >), 

3) B lU = 4> x + 0(A r2 ), B 2 u = <f> 2 + 0{h rs ~ p ), 

4) ||u - i/|| 0 = 0(h r *). 

Then, for p < s < p* , 

ll« - V llr+s + |« - V| r+a _l = 0{h r ), 

where r := min{r 1 ,r 2 ,r 3 + \ - s,r 4 }. 

4. The One-Step Scheme 

We now describe the transformation of the reduced equation (2.5) to the one-step 
scheme (2.7). We first operate on the t-th equation of the reduced system (2.5) with 

(* + -A 0 )- a ‘ + ' f (4.1) 

where A 0 is the symbol Ao(A,cu) defined in equation (2.9). Recall that 0 { < 0 and p > 0. 
Let 

i'ii ■= (<+ - A 0 )-’< +f L ii . 


16 



The resolvent equation for this new system has the additional root z — 1 + AAo with 
multiplicity \a\ + np. Since 1 + AAo > 1, no additional boundary conditions are required. 

The first p boundary conditions are modified by multiplying the fc-th boundary con- 
dition by Ao _1_Pfc and, when q is greater than p , the last q — p boundary conditions are 
multiplied by h pk ~ p+1 . Note that 

p — 1 — pk > 0, for 0 < k < p, 

p — 1 — pk < 0, for p < k < q. 

The resulting system of difference equations is elliptic of order (0, r + p) and the 
resulting boundary conditions are all of weight greater than or equal to —1. We obtain 

n 

= F/(z l/ ,w), u > 0, * = (4.2) 

3 - 1 

n 

^B'{h,y,ui)uj{x 0 ,u>) = k = l,...,n. (4.3) 

3 = 1 

This system of difference equations will now be written in a more canonical form. The 
left side of each equation of the reduced system (4.2) can be written as a sum of terms of 
the form 

P{h,y,v)6l TfcU/(*i/,w), 

where p(h,y,u>) G S r ' J+ ^ _0 ,a < ry + p and, by Assumption 2.1, 7 > 0. If 7 > 1 and 
a < Tj + p — 1 we can rewrite this term as 

p(h,y,u)6lT^~ 1 u J (x t/t ui) + hp{h,y,uj)6^ +1 T^~ 1 u[x l/ ,u). 
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In this way the order of translation 7 is reduced by one and the order of differencing 
is increased by one. Note that hp(h,y, cu) G S r ^ +P ~ 1 ~ a . By continuing in this way, we 
obtain a system where the nontrivial translation operators, i.e., 7 > 0 , operate only on 
the highest order differences. 

The resulting system may be written as 


T)6+ +P uj(x l/ , 

j=i 


«) + £ 


3= 1 


Tj+p-l 

^2 ii3a{h,y,u)8%Uj{x u ,u) 

a=0 


= £/(x„,w), 


(4.4) 


v > 0, * = 1, ..., n. 


where gij is a polynomial in T (the translation operator in the x-direction) with coefficients 
in S°, and f, ja G S T i +p ~ a . 

Before transforming the system (4.4) to a one-step scheme, it is necessary to exam- 
ine more closely the matrix of translation operators G(T) := ( gij(T )) . The extent of 

9ij{T)f>+ is at most (0, af + + p), by Assumption 2.1. Decompose G(T) as a 

sum of G\{T) and Go{T ), where Gi(T) is composed of elements g'ij(T) which, if they are 
nonzero, are of degree 7,7 = af + /?+ — <r t - — ry, and G 0 {T) has elements p°y(T) of degree 
less than 7 t y. 


Lemma 4.1 

G\{T) is equivalent through elementary row operations to a diagonal matrix, i.e., we 
can assume that 



if i = j, 


fJ-i := 


(4.5) 
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Moreover, ( 4 . 4 ) is equivalent to 

T^6^ +P ui(x u ,u) =~Y21 l, 9ijb(y,oj)T b 6^ +p Uj(x u ,uj) 

j= 1 b = 0 

( 4 . 6 ) 

n Tj+p-1 

E iij{h,y,u)6^u j (x l/ ,u)+F-\x v ,u}). 

j= 1 a=0 

1/ > 0, t = 1, n. 

By an elementary row operation we mean permutation of the rows, multiplication of 
a row by a nonvanishing function of (y,w) which is in S°, or addition of a translate of one 
row to another row. 

Proof: 

Gi (T) = (£,(/>, y ,»)rw|, 

and since det(y^) is the coefficient of the highest power of z in R(y,u),z), det(y^) does 
not vanish for any (y,w). Also, in adding to row l a translate of row m we can restrict 
ourselves to the translation operator t( q j + Jt j s then easy to see that the 

elements of Gi(T) will remain monomials in T of degree 7 ,y, although the rows, and thus 
the 7 ij, will be permuted. 

It is a standard result that by means of these elementary row operations G\{T) can 
be transformed to upper triangular form (see Gantmacher [ 5 ,p. 135 ], and since Gi(T) 
remains a matrix of monomials in T it can be transformed to a diagonal matrix. This 
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matrix has nonvanishing diagonal elements; by dividing through by the coefficients of T' T<< 
we achieve the form (4.5). 

The lemma will be proven if we can show that the elementary row operations keep 
the degree of the (t, j ) — th element of Go(T) less than 7 ij and do not create new elements 
for G\(T) from the second sum of equation (4.4). It is easy to see that by operating on 
the i-th row of (4.4) with a translation of degree (a^ — o\) — (a^“ — oy), and adding it to 
the Z-th row, the degree of the (/, j)— th polynomial in Go remains less than 7 y. The effect 
of this operation on the second sum in (4.4) results in terms of the form 

lijaT^S+Uj (x^, w), 

in the Z-th row, where 7 = (a^ - a/) - (a+ - <r t ). If 7 + a < ry + p, this term can be 
written as 

Z~»ya(l + h6 + y6*uj(x uy u), 

which, when expanded, again has the form of the terms in the second sum of (4.4). If 
7 + a > Tj + p then by repeated substitution of 1 + h6 + for T one obtains terms as in (4.4) 
with the term of highest power in T being 

T'i+ a - r >-P6? +i} u j {x l/ ,u). 

Now by inequality (2.14) a < Tj + p < at + /?t - <7, + p , so 

7 + a - Tj - p < 7 + at + fit - a - Tj < (a^ - a { ) + (/?+ - ry), 
and hence this term is also of order less than 7/y. This proves the lemma. 
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We now define the variables Wj >a (xi/, w) which are the components of the vector W in 


the one-step scheme (2.7). We set 


Wj, a [x v ,u) := Al i+P 1 °^uy(s„,w) for a = 0, ...,Ty + p - 1, 
If fij > 0, we put 


Wy >a (x„,a;) := T a (Tj+/i a) 6j +p 1 Uj{x v ,w), for a = ry + p, ...,ry + p - 1 + /zy. 

Thus there are a total of |/z| + |r| + np components in W. 

The equations comprising the one-step scheme are: first, if ry + p — 1 > 0, 


Wj-, a (x u+u u) = Wj t a[x v ,u) + hA 0 W j>a+ i(x vt u) for a = 0, ...,ry + p - 2, (4.7) 

and, if /zy > 0, 

Wj t a{x„ + I,w) = Wj t a+l{x v ,u) t for O = Tj + p - 1, ...,Ty + p + flj - 2. (4.8) 

For a = Tj + p + Hj — 1 we use equation (4.6) to obtain an expression for Wj fa {x„,u). 
Multiplying equation (4.6) by h we have 




T^S^-'u^.u) -it, II - T^-'u^u)] 

3 = 1 b = 0 
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n T j + P —2 

+E E ™ijbWj t b[x vy u) + hF"(x Ui w), 

j= 1 6=0 


where rhijb — ~hl”j b A 0 * :+p 1 ^ € S°. Replacing T b 6^ +P 1 uj(x u , u>) by Wj^+rj+p-i in 
the above, we have (with a := r t - + p - 1 + /x t ) 


w) — (4.9) 

n Mj-1 

“ X^ X^ + - ^J,r,+p-l+&(z l /,w)] 

3 = 1 6=0 

n ry+p— 2 

+E E rhijbW j, b {xi,,u) + hF"(x y ,w). 

j = 1 6=0 

The equations (4.7), (4.8), and (4.9) together give the one-step scheme (2.7). We also 
will write equation (2.7) as 


Wi, a {x v+U u3) = M ij*b{y,u)W jtb (x u ,u ) + h7i[x uy u), (4.10) 

3 , 6 


* — 1 > •••> n > a 0 , ..., T( + P — 1 + 

Note that the matrix iVf(t/,w) as constructed is independent of h, depending only on 
y and w. This follows from equation (2.4) and the means used to define the W iia (x u ,u). 

Note that hAo{h,u) = \Ae£ i sin 2 (|w t ) is independent of h, depending only on w. 

The boundary conditions (4.3) can be written in terms of the variables W ifa as 


Bi{y,u)W(x 0 ,u) = $i(w) 
B 2 {y,u})W(x 0 ,u>) = 


(4.11) 
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where Si and 82 consist of the first p and last q — p boundary operators, respectively, as 
follows. Each of the terms in the sum (2.6), describing the boundary conditions for the 
reduced problem, can be written as 


Tj+Pk 

BkjUj = ^2 b k,j,a{ fl ’y’V)ho +Pk ~' a 6+Uj, 

a = 0 

where the 6fc,y, a are symbols of order 0 in uj. The first p boundary conditions are scaled 
by Ag -1-P \ as mentioned earlier in this section, and thus 


Tj+Pk Tj+Pk 

K~'~ n B km = Y 

a=0 a = 0 

This defines the elements of Bi(y,u;). The components of $1 (w) are A£ 1 pk <f>k{w). 

The elements of B^y, w) are obtained by scaling with h ph ~ p+1 . Thus each term in the 
sum in (2.6) for k from p + 1 to q becomes 


Tj+Pk 

/.«-*>+■ Y b >.kA 


’j+Pk-Ora . 
0 °+ u J 


Tj+p-\ 


a= 0 


a = 0 


T j+Pk 

a=Tj+p 
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The first summation can be expressed as 


Tj + P- 1 

E b i.kA^oY k -^'Wi, c 

a= 0 


( 4 . 12 ) 


and the second as 


Pk—p+i 


" a / \ 

E t>iX' + r i+f -,(hA 0 ) n -' + ‘-‘ E (“ 

a=l c =0 


This second summation is further expressed as the sum of several terms, 


pfc-p+i 


min ) 


X* ^7,M+^+p-l(^Ao) P * P+ X ( _ ) ( l) a C ^ / y,c+r i -fp-li (4-13) 


a=l 


c =0 


and 


PJfc-p+l a 

^i»fc»o+r,-+p-l(hAo) W ^ +1 a 
“I" 1 c=/xy + l 

where this last expression is taken to be zero if pj is greater than Pk~P+ 1 . The expressions 
T c ^kFy.rj +p-l+Pj are replaced by repeated use of the equation TW = MW + hi. This 

defines the elements of 82(1/, u>) and the components of $2^) are h Pk ~ p+1 plus terms 

from hi. 

We now state and prove two lemmas about the matrix M. 

Lemma 4.2 



(“) (-l)-T (4.14) 


det ( — — — ) = h M\^ +np Hh, y,w)c(y,u), 
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where L(A,t/,w) = det(lij(h,y,$)) for f = (0,u;), c(y,u>) / 0 for all u, and c(y,u>) is 
independent of h. 

Lemma 4.3 


iit^m < ca^, 

for some constant C. 

Proof of Lemma 4.2: 

There are three types of rows in the matrix (M — I) jh. Thinking of W as a doubly 
subscripted vector, the first type is obtained from equation (4.7) and contributes one 
nonzero element, Ao, in row (i, a) and column (i, a+l) for a = 0, ..., T t +p— 2, 1 < % < n. The 
second type, obtained from equation (4.8), contributes two nonzero elements —A -1 and A -1 
in row (t,a) and columns (i,a) and (t,a+l) for a = r,-+p— l,...,r< + p+p,— 2, 1 < t < n. 

We now evaluate the determinant. First, each row corresponding to an equation of 
the type (4.7) has only one element, Ao, in the position ((j,o),(j,a + 1)),1 < j < ra,a = 
0, ..., Tj + p — 2. Evaluating about these rows gives A^ +n ^ -1 ^ times the determinant of 
the remaining rows and columns. The second class of rows corresponds to equation (4.8). 
These have the elements A -1 in the position ((j,a), (i, a)) and —A -1 in ((j,a), (j, a + 1)), 
1 < j < n,a = Tj + p — 1, ..., Tj + p — 1 + fij — 1. However, in evaluating about the first class 
of rows the column (j, Tj + p— l) is eliminated, so by evaluating the determinant remaining 
after the first reduction there results a factor A -1 times m,y 0 (w) in the (t, r, + p — 1 + p,) 
row and the (y,0) column which is, by definition, equal to / t 'yo(A, w)Aq Thus, from 

the definition of the /,-yo resulting from the operators (4.1), the determinant of (M — I)/h 
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IS 


±h l^lA| ) a|+np det(/^(fc, y,u>))/det(^ y (w)), 

where g[j(w) are the coefficients of Gi(T). Since det(yU(u>)) does not vanish for any to, 
this proves the lemma. 

Proof of Lemma 4.3: 

Using Cramer’s rule to evaluate the inverse matrix of (M — /)/ h one sees that each 
minor is bounded by either Aq 1 or h times the determinant of (Af - I)/h. Since h<c/ A 0 
the estimate follows easily. 

Lemma 4.4 

For to — 0, 1 is a semi-simple eigenvalue of M. That is, there are exactly |r| + np 

linearly independent vectors V j ’ b ,j = l,...,n,6 = 0 ,...,ry + p - 1 such that MV*’ b = V>' b . 
Moreover, these eigenvectors are given by 


yJ> b 

i,a 


f>i,jt>a,b, if 0 < b < Tj + p - 2; 

+ if o, > T{ + P — 1, 6 = Tj + p — 1; 

0, otherwise. 


Proof: 

We consider the equation MV = V. From the (4.7), we see that (MV) iia = V ij0 for 
0<a<Tj + p-2 and from (4.8) (MF) t>a = U tja+1 for ry + p - 1 < a < ry + p + /xy - 2. 
Thus if V is an eigenvector of M with eigenvalue 1 then V, i0 = for these 

components. In (4.9) note that A -1 m t) y ( & £ S 1 and hence vanishes at w = 0. Thus for an 
eigenvector, (4.9) reduces to {MV) i>a = V i>a for a = + p - 1 + fr. It is now easy to check 
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that the |r| + np independent vectors V 3,b given in the lemma are indeed eigenvectors. 
Moreover, by Theorem 5.3 the multiplicity of the eigenvalue 1 is |r| + np, hence these are 
all the eigenvectors and 1 is a semi-simple eigenvalue. 

Lemma 4.5 

If MW = W for w = 0, then B 2 W = 0. 

Proof: 

If u> = 0 then Ao = 0 and all the coefficients in (4.12) vanish since pk > p for k > p 
by Assumption 2.2. Since MW = W the sum of (4.13) and (4.14) gives the coefficient of 

b],k,a+Tj+p - 1 3 s 


min(a,^y) 


nyu’iH’j / / \ a / \ 

E ! E w<-i) CW 3, T j+p-\+nr (4*15) 

c=0 W c= N + \ W 


Again since MW = W the equation (4.8) implies that 


Wj^ a — Wy,a+ 1 > for O, — Tj + p 1, + p + Pj 2, 


and hence (4.15) is equal to the quantity 



7 


which vanishes identically. This shows that B 2 W = 0, as asserted. 
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5. Construction of H(y,u) 


We now construct the matrix H(y,w) to satisfy inequalities (2.9) and (2.10). 
Theorem 5.1 

There exists a bounded Hermitian matrix symbol H(y, u) £ S° which is a C°° function 
of ( y,u ;) for y £ R d_1 and u £ R i_1 \ {0} such that 


a) M*HM - H > c 0 h\ 0 

b) W*HW > Cl {n\W + \ 2 - |W_| 2 ), 


(5.1) 


where W + = W+(w) is the projection ofW into the span of the generalized eigenvectors of 
M whose eigenvalues are of modulus greater than unity, and W- = W_(u>) is the projection 
into the span of the generalized eigenvectors of M whose eigenvalues are of modulus Jess 
than unity. The constants co,ci and tj are all positive and, moreover, tj can be chosen 
arbitrarily large with c 0 ,c i, and H depending on jj. 

The proof of Theorem 5.1 will be postponed until after we prove Theorems 5.2-5.4. 
We begin by constructing H in a neighborhood of iv 0 ^ 0. 

Theorem 5.2 

There exists a smooth, bounded matrix function P(oj) defined in a neighborhood U 
of u > o ^ 0 such that 


P(w)M(w)P(w) 


i = ( M+{u) O \ 

^ O M_(u>) ) ’ 


(5.2) 
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where M + (oj) is a m+ x m + matrix and M-(u>) is a m_ x m_ matrix. Moreover, for u> EU 


M+(u)M + (w) > 1 + c 0 , 
and AfI(w)M_(w) < 1 - c 0 . 


(5.3) 


m+ is the number of eigenvalues of M(oj) of modulus greater than unity and m_ is the 
number with modulus less than unity. 

Proof: 

By lemma 4.2, M(u>) has no eigenvalues of modulus one for u> ^ 0. Therefore it 
follows by standard linear algebra that M(w) can be transformed to the form (5.2). That 
the inequalities (5.3) can be satisfied is also a standard result, see Gustafsson et al. [6]. 


For u)q ± 0, H(uj) can be constructed in the neighborhood U of wo as 

where 7 m+ and I m _ are the identity matrices of order m + and m_, respectively. It is 
easily checked that H(u) satisfies (5.1) in U. 

We now consider M(u>) in a neighborhood of a>o = 0. 

Theorem 5.3 

For |cc7 1 near zero the eigenvalues of M (cj) separate into three distinct classes. These 

are: 

1) There are |o| + np eigenvalues with k = 1 + HAq. 

2) There are |r| — |cr| = 2p eigenvalues which are the roots of the equation 

det(lij(h,$)) = 0, 
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with 


f = (Ci,w) and k = exp (»fi), 

where 

k . = 1 ± 0(/iAo). 

Moreover, p of these eigenvalues satisfy 

|/c| < 1 — c/iAq, 


and p of them satisfy 

|/c| > 1 + c/iA 0 , 

for some positive constant c. 

3) There are |/x| eigenvalues which satisfy 


|«| ~ 1 | > 6 > 0 . 


Proof: 

The proof depends on the equivalence of the one-step scheme to the original reduced 
equations (2.5). The first class of \o\ +np eigenvalues with k = 1 + hA 0 are due to operating 
on the reduced system with ( S + - A 0 ) _£T *' + ^ (see 4.1). There are a total of 

n 

+ p) = M + n p, 

i=l 

such eigenvalues. The third class of eigenvalues is determined as follows. At wo = 0 the 
eigenvalues of M are the roots of 
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det( 9i ,'M)(*-l)' r|+,W ’ = 0, 


by equation (4.4). Notice that the / tJa in equation (4.4) vanish at u 0 = 0 since they are 
symbols in S rj+ ^ -1_a with only highest order terms. Since the system (4.4) is elliptic of 
order (0, r + p ) the root k = 1 cannot have multiplicity greater than |r| + np. These are 
the roots of class (1) and (2). Thus the roots at uo = 0 which are distinct from 1 are all 
the roots of 


det (#,-(«)) = 0, 

of which there are \p\. By Lemma 2.1, |k| is not equal to 1. By continuity there exists 
some neighborhood of u>o = 0 and constant for which the inequality in part (3) holds. 

The remaining |r| — [<r| eigenvalues are easily seen to be related to the associated 
system of partial differential equations (2.15) . If we define k(6) by 

det(/ f y(y,«,0)) =0, 


then 

= 1 + M*(o) + °(M 2 )- 

M 

By the supplementary condition (Assumption 2.3) we easily obtain the inequalities in 
part (2) for the p eigenvalues less than and greater than one in modulus. 

Theorem 5.4 

For u in a neighborhood of zero there is a continuous nonsingular matrix Q(co) such 
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that ||Q(w)|| and \\Q 1 (w)|| are bounded and 

L 0 (u ) O O O \ 

o N °0 ] L°w) O - ( 5 - 4 ) 

Vo O O JVjH/ 

where Lq,L\,Nq, and N\ are square matrices. Moreover, 

fl) L 0 io ^ 1 — <5, 

6 ) NqN 0 >1 + 6 , 

(5.5) 

c) ZrjZ /1 <1 — Cl/lAo, 
rf) iVjiNTj > 1 + cxAAo, 

with Z/ x (0) and Ni( 0) being identity matrices. The dimension of Li is i(|r| - |a|) = p and 
that of Ni is |(|r| + |a|) + np. 

Proof of Theorem 5.4: 

For w in a neighborhood of zero the eigenvalues of M(u) separate into three classes: 
those which are strictly less than one in modulus, those which are strictly greater than one 
in modulus, and those which are not bounded away from one in modulus. Thus Q(u) can 
be constructed so that 

( L 0 O O \ 

Q{u))M{u))Q- 1 {uj) = O N 0 O , 

\0 O Mi J 

where the three matrices on the diagonal correspond to these three classes of eigenvalues. 
Moreover, Q(w) can be constructed so that L 0 and N 0 satisfy the inequalities of the 
theorem. 
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By Schur’s theorem we can assume that Mi(u;) is in upper triangular form with the 
first |(|r| — |cr|) diagonal elements being those eigenvalues which are less than one in 
modulus for u> nonzero. There are |(|r| — |a|) = p such eigenvalues by Theorem 5.3. Thus 



M 2 (u>)\ 

n 2 H ) 


We now show that there is a matrix D(w) such that 


[I d\(l 2 m 2 \ ( I -d\ = (l 2 o\ 
\0 I ) \ o n 2 )\o I ) \0 N 2 )' 


with ||D(w)|| bounded for u near zero. 

The matrix D(u>) is the solution to 

L 2 (u)D(oj) — D(oj)N 2 (u) = M 2 (w). 


In order to prove the existence of D(u>) we give the following 
Lemma 5.1 

There is a unique solution to 


L 2 D — DN 2 = M 2 , 


with 

IPII < jr||Mj||/fcA„, 
where K is independent of h and u>. 
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Proof: 


By Schur’s theorem there are orthogonal matrices Ox and 0 2 such that L := Oi*L 2 Oi 
and N := O 2 N 2 O 2 are lower and upper triangular, respectively. With M := OlM 2 0 2 and 
D := 0\D0 2 , the equation becomes 


{La - N kk )D ik = M ik - Lijbjk + £ D t: N jk . 

j<i i<K 

This is a recursive formula for the elements of D in the order D n , D 12 , b 2 i, D a b- 
By Theorem 5.3 (La ~ N kk ) = O(hAo) and the estimate follows easily. This proves the 
lemma. 

It remains to show that ||A^ 2 [[ < cKAq. We will need the following lemma, whose proof 
will be given after the completion of the proof of Theorem 5.4. 

Lemma 5.2 

Let x(e) be an upper triangular matrix depending on e for 0 < e < e 0 . If 

II < c/e , 

and 

c\e < |x,-,-| < C\e for 0 < e < ep, 

then 

l*»y(e)| < c £, 

for some constant C. 
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Now with e = hA.o and x{e) = M\{w) — I it follows from Lemmas 4.3 and 4.4 that 
HM 2 II < c/iAo . Hence D(w) is bounded in norm for u> near zero. 

An argument similar to the above shows that the off-diagonal elements of L 2 (u>) and 
N 2 (u) are all of order hAo. Thus there are constant diagonal matrices L>i(w) and D 2 [u) 
such that for 


L\{w) DiL 2 (u>)D 1 x , 

JVi(w) := D 2 N 2 {u)D^, 

we have 

L\L\ <1 — c/iAo, 

N{N\ < 1 + chA 0 . 

This completes the proof of Theorem 5.4. 

Proof of Lemma 5.2: 

The proof is by induction on j. Suppose (x,y) < C(l)e for j = i + /, 0 < / < m — 1 . 
For j = i, \xij\ < Ce by hypothesis. Let z,y := (x _1 (e)),y, then 




-1 


i+m-1 


y =l - 


( x ij z j,i+m)i 
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and hence 


t'+m-l 

|^*,t+m| < C\£ (C(j)e)(c/e) < C{m)e, 

j=* 


which proves the lemma for C equal to the maximum of the C(l). 


For M(w) in the form (5.4) satsfying (5.5), we construct H as 


/-I O O 0\ 

O rjl O O 
OO-IO' 

Vo o o r)ij 


and then 

H{oj) := Q’(u)HQ{ w). 


It follows that 

M*HM - H = Q*(M*HM - H)Q 
> cohAoQ* Q 
^ c 0 /iA 0 . 

This satisfies (5.1a) locally; similarly (5.1b) is satisfied. 

We have constructed H(y,u>) satisfying (5.1) in a neighborhood of each u. Since the 
set of u with |w,| < tt is compact we can choose a finite set of these neighborhoods which 
cover R d 1 , and through the use of a partition of unity we can construct H (u>) for all cv. 
This proves Theorem 5.1. 
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We now show that the relation (2.10) is satisfied for some positive constants c\ and 
C 2 . By equation (5.1) we need only show that 

|BW_| 2 > c\W_\ 2 , 

for some positive constant c where is defined as in Theorem 5.2. For then 

W*HW + c x |BW| 2 > c x ri\W+\ - Cl |BP^ + | 2 + (c x c - c x )|W_| 2 

> c 2 (\W + \ 2 + \W_\ 2 ) 

> c 2 \W \ 2 , 

for Ci and r] chosen large enough. Thus we need to prove 
Theorem 5.5 

If the Complementing Condition is satisfied then there exists a positive constant c 
such that 

|B(w)W_| 2 > c|W_| 2 , (5.6) 

for all <jj G R d-1 f and all W- defined as in Theorem 5.1. 

Proof of Theorem 5.5: 

Suppose for some wq with u>q ^ 0 

5(w o )W_ = 0 for \W_\ = 1, (5.7) 

where W_ is in the generalized eigenspace of M (w) whose eigenvalues are of modulus less 
than unity. Then by setting 

W 0 := W- 
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w v+l := M(uo)W l/ 


v > 0, 


W w is an eigensolution of type 1. W u tends to zero as u increases since W_ is in the 
generalized eigenspace whose eigenvalues of modulus less than unity. Thus (5.7) violates 
the Complementing Condition. 

Next suppose that (5.7) holds with u>o = 0. It follows that 

B 2 (0)W- = 0, 

and we see that W u , generated by M(0) as above, corresponds to an eigensolution of type 
III. Therefore (5.7) violates the Complementing Condition when w 0 = 0. Now consider the 
possibility of the existence of sequences and such that 

—*■ 0 as k — * oo, \W±\ = 1, 

and 

> 0 as k — ► oo. 

We may assume that there exists a such that 

and u k l\u k \->e as k - oo. 

Set W® := Wjj + Wjjj where Wjj is in the linear space spanned by the eigenvectors of 
M with eigenvalue 1 and Wjjj is in the span of the eigenvectors with eigenvalues strictly 
less than 1 in magnitude. By lemma 4.5, B 2 W^ = 0, and, since S 2 W° = B 2 W^ n = 0, the 
nonexistence of eigensolutions of type III implies Wfjj is zero. Therefore, without loss of 
generality, we can assume 
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w! = ittfaW’ 

6=1 

where the /Cj,(u;y) are the eigenvalues of M(u>y) with /Cf,(wy) — »• 1 as wy — * 0, / fc fc — ♦ /°, and 
|/C{,(u;y)| < 1 for wy 7 ^ 0 . Define the vector-valued function Z(x) by 


Z(x) ■- £/?«->** 

6=1 


where 


rjb := lim 

k— MX 


1 - /Cfc(cjy) 


IWv 


By the structure of the eigenvalues of M at cu = 0 (Lemma 4.4), we have that Z Jta (x) = 
Zj, T ,+fi-i[x) for Tj + p — 1 < a < Tj + p — 1 + /tty. From (4.7), for 0 < a < Tj + p — 1 we 
have 


Wj, a,v+\ Wy, a ,i/ — hA.oWj )a +i )V , 


or 

p 

fb,j,a K b{ K b ~ l) = hAoVFy^+i^. 

6=0 

Dividing by |wy| and taking the limit we have 




6=0 


6=0 


hence 



'Zy.a+iCz) 
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for 0 < a < Tj + p - 1. Thus 


Q 

z iA x ) = (^) az j.o( x ) for ° = °> -.»>• + P - 1. 

Let Zy(x) be Zy,o(x). Finally, from (4.9), we obtain 



j = 1 6=0 

By reversing the transformations done after equation (4.4) (see Lemma 4.1) we have that 
this system is equivalent to one analogous to (4.4), namely, 


n r j+p— 1 


EMo,M(^r’ + *2,M+E E 


y=i 


y=l a=0 


By eliminating the factors corresponding to (4.1) we obtain the associated 

system of differential equations (2.17) and the vector function Z(x) = (Zy )0 (x)) is a solution 
to (2.17). 

By a similar reduction one can show that Z(x) satisfies 8iZ( 0) = 0. The nonexistence 
of eigensolutions of type II implies that Z(x), hence Z(x ), is zero. Thus we have, for 
|u>_| = 1, that \B(uj)W-\ > c which proves Theorem 5.5. 
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6. The Regularity Estimate 

In this section we prove the regularity estimate (3.3) using the inequalities (2.9) and 
(2.10) and Garding’s inequality (Bube and Strikwerda [3]). We state the final estimate as: 

Theorem 6.1 

If the elliptic difference scheme (2.1) with boundary conditions (2.2) satisfies the 
Complementing Condition and Assumptions 2.1, 2.2, and 2.3, then the following regularity 
estimate holds for p < s < p* and h sufficiently small, 


imi?+. + m? + ._ 4 < + 1 + ii^ii a-. + Huns). (6.i) 

where t = p + | [2(s - p )\ . 

The first step to proving Theorem 6.1 is to prove estimates on the tangential differ- 
ences. We first define the norms 

n«ii?+<,r+. = ||A , ' + *— *>>||?, 

j— 0 a=0 


and 


u 


lr+i,r+r 

j —0 a =0 


which limit the normal differences, i.e. those with repect to x, which are included in the 
norm. 
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Our first result is 


Theorem 6.2 

If the elliptic difference scheme (2.1) with boundary conditions (2.2) satisfies the 
Complementing Condition and Assumptions 2.1, 2.2, and 2.3, then the following estimate 
holds for any real s. 


+ |fe p - f+ V 2 lL ji+ ; + 

+ ll“ll?+p-l,0 + Mj+p-i.-l)' 


( 6 . 2 ) 


Proof: 

We have that 

OO 

ii/iis, 

u=0 

where \fi,\o,s is the Sobolev norm of order s on the tangential variables only at x v . Using 
equation (2.9) and G&rding’s inequality on the tangential variables, we obtain 


IIW'llo,. < C YJXM'BM - 


u=0 


where the tilde indicates the pseudo-difference operator corresponding to the given symbol. 
For solutions of (2.7) this gives 


Iloilo, /< C(-QVo,HWo) 0 ,.-t+e 


h + -mile-! 
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Equation (2.10) and G&rding’s inequality again imply 


+ IH'olJ,..! < C(\BW 0 \l . + ellWII?,, + - ||71lo, + IIM'llo,,-* + l^ol?,-,). 


Taking e sufficiently small we obtain 


+ IH'olS,,.! < + Iloilo, «-i + IIH'llS.-i + 


Also, for er > 0 there is a c s>e such that 


2 

0,s — 


1 

2 


< c\m\ 


2 

0,3 


+ c.,,\\W\\ 


2 

0 , 0 > 


and similarly for the boundary norm, e.g. Thomee and Westergren [ll], or Bube and 
Strikwerda [3]. Thus we obtain 


IIH'lliv + IH'oi;,.-* < C(|BH'ol?,,_i + imiS,.-, + IIW'llg.0 + (6-3) 


The estimate (6.2) for the original dependent variables is obtained from (6.3) by replacing 
the Wj t a by the equivalent expression in terms of the uy. The norms and scaling of the 
boundary data <j>\ and fa in (6.2) are the result of the modifications to obtain the one-step 
scheme as described at the beginning of Section 4. 

Estimates on the Normal Differences 

We now consider estimates on the normal finite differences, and begin by stating 
interpolation results for the normal differences. 
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Lemma 6.1 


For integers t, s with t < s and e > 0, there exists a constant C £ such that for any 
discrete function i ’(x r ,y M ) 


\ 6 . 


+ v lo,r 


< 4K +1 



+ C e \\6lv\ 


2 

0,r 4- 


1 

2 


and 


||£ft>||o,r < vjjo ,* — 1 + C e ||<5^v||o, r +l- 


Proof: 

By summation by parts, 


OO 

(A r fJi>(x 0> w)) J = -h\' Y, w)) 2 

u=0 

oo 

= ~ h ^( A " (r -" ) <5; +1 t;(x l/ , W )) A i( r +^ ) (6; V (x l/+ i,u,) + 6^v(x Vt u)). 

u—0 


Hence 


I s + v 


2 

0,r 


< 4CWI 



+ C,\\S‘ + v II 


2 

0,r+§- 


Then by using an interpolation inequality (e.g Thomee and Westergren [ll]) the first 
inequality of the lemma is easily proven. The second estimate is proved in a similar 
manner; the proof will be omitted. 

For the case n = 0 the restriction on the differences in x in Theorem 6.2 is removed in 
precisely the same way as it is for the partial differential equation. From equation (4.6) we 
have that 6 + +P «y(ii/, u) is equal to a sum of lower order differences of uy and F"{*u,u). 
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We easily obtain 


||£^ +P «||o,s— T- p+1 < C’(||«||r+p-l,r+s+l + ||F||. + ,-.). (6.4) 

The use of the interpolation estimates in Lemma 6.1 can then be used on the norms of u 
on the right-hand side of (6.2) to give (6.1) for s - p in the case that n = 0, which implies 
that q = p. By operating on equation (4.6) with 6+~ p and using an estimate similar to 
(6.4) along with Lemma 6.1 we obtain (6.1). 

The case \i 0 is more difficult because one must solve for 6^ +p Uj(x l/ ,u) from 

equation (4.6). We now show how this is to be done. It essentially involves solving an 
elliptic system of order (0,0) for the 6+ +p Uj(x u ,u). 

Consider the system of equations (4.4) and write it as 

G(u,T)6+ +p u(x l/ ,u) + K(u)u(x„,u>) = F'(x v ,w), 

where 6+ +p u(x u ,u) denotes the vector with components 6jf +p Uj(x u ,u), j = l,...,n, and 
G{u,T) is the matrix of translation operators <7,y(u>, T). The term K{w)u{x v , w) contains 
all the differences of the Uj(x t/i u) with respect to x of order less than ry + p. The operator 
G is an elliptic operator of type (0,0) by the Resolvent Condition. 

As in the proof of Lemma 4.1 the matrix G(u, T) is upper triangular and without loss 
of generality we can assume that the rows are ordered so that the degree of < 7 , )t (u>, T) is 
greater than the degree of &y,y(w,T) if i < j, i.e. /«,- > f/y if i < j. Let n' be the integer 
such that m = 0 for i greater than n', with n' = n if // n > 0. For i > n', <5^u t - can be 
expressed in terms of lower order differences of the uy as in the case h — 0. 
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As in section 4 we construct a one— step scheme for W the vector whose components 
are T 6+ ~u J (x l/ ,u>), for 0 < / < pj — 1 and 1 < j < n' . This scheme can be written as 

W«/+i = MW u + g u (6.5) 

where contains both the function F , and difference of uj with respect to x of order less 
than Tj + p. The boundary conditions for this scheme can be written as 

CW o = $ ( 6 . 6 ) 

where contains both the data $2 from (4.11) and differences in uj(xo,oj) of order less 
than r + p, and C has been obtained from (3.2) by scaling with h p ~ p analogous to the 
scaling at the beginning of section 4. 

The matrix At is of order |//| and it is easily shown that for u near 0 its eigenvalues 

are those of class 3 as given in Theorem 5.3. Thus we can construct a symmetrizer ){ such 

that 

H > c 0 (6.7) 

as in section 5. By the nonexistence of eigenvalues of type III, there is a neighborhood for 
u i near 0, such that 

)i J rc l C*C>c 2 (6.8) 

for some positive constants Cl and c 2 . Note that if q = p then C is taken to be zero, and 
^ is a positive definite matrix. The relation q = p results from all the eigenvalues of .M 
being larger than 1 in magnitude. The inequalities (6.7) and (6.8) are analogous to (2.9) 
and (2.10). 


46 



We now extend the estimate (6.2) to include normal differences of r + p. Since (6.8) 
only holds for small values of u we employ a cut-off function ip(u>) with i>(u) = 0 for 6q < 
Moo < 7r, xjj(u) = 1 f° r M°° < 6o/2 and extend xp(u) periodically for all u e R d_1 . Let 
"Vv be xl){uj)M> u . Analogous to the proof of Theorem 6.2, but using the Garding inequality 
proved by Lax and Nirenberg [7], we have 

OO 

imio.r < Ch £{((£*#/« - X)Wv)o,r + ch{V vt V v ) o.r}, 

u-0 


< c{—h{v 0 , ^/Vo)o,r + e ll^llo,r + ^e|l5llo,r + c Ml^llo,r}- (®*9) 

By (6-8) and the definition of C and Q we have for small h 

PUS,, + fe|Vo||, r < C{h\CVo\lr + H5llo,r} (6-10) 

< C{h\u\ 2 r+p _ Ur+l + |k'-' + 4 2 |? + ||F||? +m + |M|? +? -, 

Since hA is bounded, we have 

h*\u\ T +p-l,r+l < C|tt| r+ ^_l, r+ i. 


For Moo > ^o/2 we use equation (2.7) as follows. We have 

i(W(x„ + i,w) - W(x l/ ,u;)) = - I)W(x U ,u) + 7 {x w ,uj), 

and estimating only the differences of the components W t)TV+ p_i of W we obtain 

\6+ +p u(x„,u)\ 2 < ^|W(x t/ ,a;)| 2 + |7(*„,w)| 2 < CA%\W[x Vi u)\ 2 + |^(*«/>M| 2 > ( 6 - n ) 
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since h 1 < C(S 0 )A 0 for Iw^ > 6 0 /2. 


Combining the estimates (6.10) and (6.11) we obtain 


E < E II* ♦ *? + 'Wll? + E 11(1 - *) * 

j =1 j'=i y=i 

< c{\h*-^<h\i + + n«ii; + ,. 1>T+Wr + fci«i; + ,.,, 


Then, using estimate (6.2) and Lemma 6.1, we obtain 


ll U llr+p,r+ a + |«|? + p_ 1(r+a _i <C{\4>i\]_ p _ h + \h p ^ + *<f>2\l-p+\\F\\g_ a . + ||u||o}, (6.12) 


which proves (6.1) for p < s < p + A. 

To prove (6.1) for larger values of s we we obtain a sharper estimate than (6.10) for 
the boundary terms, i.e. |Vo|o,r> as follows. Since the eigenvalues of X are bounded away 
from 1 independent of h, there is a number (3 with 0 < /? < 1 such that for any eigenvalue 
A of X 

c<\m\ 2 -i\, 

analogous to the estimate for the third class of eigenvalues of M in Theorem 5.3, and 
moreover, |/?A| > 1 if |A| > 1. Thus the matrix 1/ can be constructed to satisfy 


c 0 <PM*HMP-)t 

analogous to (6.7). The variables satisfy 

0 v+1 Wv+ 1 = {0M)p u iv v + p u+1 q u 
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by (6.5). Similar to (6.10) we have 


U0V||§, r + A|V 0 |2, < C{A|CVo|g, r + II/J5II?,,}. 


(6.13) 


and we estimate the last term as follows. We have 


OO 


W9\\l,r = l'E^\S\lr 


«/=0 


h 

l-p 2 


u = 0 


h 2 

1-/3 2 


OO 


s . 


2 

0 ,r 


h 2 

1-P 2 


(\Bo\ 2 +l3 2 J2^ S +M,r) 

u = 0 


< 


h 2 /3 2 

1-P 2 


OO 


^ ' l^+l^.'lo.r 


< Cfc||5llo,r+i IK+Sllo,,- 


1- 

2 


Thus we have from (6.13) 


l^olo.r ^ ^(l^^olo.r + i + ^11^ llo.r+l ) 


and thus, with r = s — 


u 


T + p,T + S 


_i < C(|u 


T+p— l,r+s — § + ^IMIr+p.r+a "l* ^IMIr+p— l,r+a 


+ \h P P 4>2? s _p_k + |^llj_ p _i + ll-f’lia-a) 


This inequality with (6.10) and (6.12) gives 


U \\ 2 T + P,T+S + M? + p >r+s _l < C{\<j) X \ 2 a _ p _y + I h P P 4>2\]_ p _L 


+ M2-<r+IMlS). 


(6.14) 
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which proves (6.1) for p < s < p + 1. 

Estimates for higher normal differences are obtained in the following way. The equa- 
tion for 6+Wv is 

6 r + W u+1 = M6 r + W, + 6 r + g u , (6.15) 

where for simplicity we have assumed that M is independent of x. If M is not independent 
of x then the right-hand side of (6.15) would contain lower order differences of W, which 
would not effect the final estimate, but would complicate the following formulas. 

For integer values of r with 0 < r < p* - p we have boundary conditions 


c (r) <5;w 0 = £ (r) 


(6.16) 


analogous to (6.6) where now C (r) is scaled by h p ~ p ~ r . In the same way that equations 
(6.12) and (6.14) were obtained we obtain 


Mlr+p+r.r+a + I u lr+p+r- l )T +a_ £ — C(l^l| e-p-i + \h P P 2 <f>2\\-p- r 

+ IMU + IMI2), 


(6.17) 


and 


MI?+,+,,,+. + \ u \r+p+r,T+s—i < C( |^l|L,-J + 1 


+ ll-f’lls-CT + ||U||S). 


(6.18) 


These estimates, with r + p<s<r + p- f- ^ and r + p + ^ <s<r + /j + l, respectively, 
prove Theorem 6.1. 

We now consider estimates on differences of order p* or higher in the case that q > p. 
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Theorem 6.3 


If the elliptic difference scheme (2.1) with boundary conditions (2.2) satisfies the 
Complementing Condition and Assumptions 2.1, 2.2, and 2.3, then the following estimate 
holds for s > p* 


Ml? + . + 2 + r—HM + l|f 111-. + IMIS)- (6.19) 


Proof: 

As in equation (6.9) we obtain, for h sufficiently small 


ii*;wns, r < c{-m«;w 0 ,«;wo) + ii«;$u§,,}. 


Let p' = p* — p — 1, then we then use equation (6.5) to obtain 


(6.20) 


, r-p'-l 

S r + Wo={h- 1 {M-I)) r ~ P 6iWo + h-^-^ £ Pv{M)6”g 0 


where p u (M) is a polynomial in At. We then have, using Lemma 6.1 on Q, 


r—p' — l 

fc|«;Wo| 2 <CA- 2 <'-' / -i)(|«''Wo| 2 + E K£o| 2 ) 

u=0 

<ci>- 2 ('-''-i)(|<jw„| 2 + ||5||L f . ), 

< Ch- 2(r - p '~ i) (|u|l +f ._i + |M|l +( ,--l,, +( ,-+r). 


from which the estimate (6.19) follows easily from (6.20) with s = r + p. 
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To complete the proof of Theorem 3.1, we state 
Theorem 6.4 

If u(x l/ ,y fl ) is a solution to the system (2.1) with boundary conditions (2.2) and 
Assumptions 2.1, 2.2, and 2.3, are satisfied, then the regularity estimate (6.1) holds for 
each s, with p < s < p* and h sufficiently small, only if the Complementing Condition 
holds. 

Proof: 

Assume the Complementing Condition does not hold. If there is an eigensolution of 
type I or type III, then as in the proof of Theorem 5.5 we can construct a solution of 
(2.1) and (2.2) with homogeneous data for any h > 0. Since the eigenvalues, /C{,(u;), are 
bounded away from 1, the norms ||u|| r+a and |u| r+s _i will be 0(h~ 7 8 ) as h tends to zero 
even though ||u|| 0 is O(l). Thus the estimates (6.1) and (3.3) can not hold. 

If there is an eigensolution of type II, then the regularity estimate analogous to (3.3) 
fails to hold for the solutions of the associated system of differential equations, (Agmon 
et al. [l]). The eigensolution is a solution to the difference equations with inhomogeneous 
data, i.e. the “truncation error”, which tends to zero as h tends to zero. It is easily seen 
that the estimate (3.3) can not hold for sufficiently small h. This proves Theorem 6.4. 


7. General Domains and Lower Order Terms 

We now discuss the modifications required to handle the cases when we have domains 
which are not a half-space or the equations have lower order terms. For a bounded domain 
fl with a smooth boundary we assume the grid is boundary-fitted. This means that for each 
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boundary point there is a neighborhood which can be smoothly mapped onto a portion of a 
half-space with the grid being mapped onto an orthogonal grid on the half-space. Through 
the use of such mappings and a partition of unity one can obtain regularity estimates up 
to the boundary on the boundary-fitted grid on Q. 

Lower order terms cause no problems unless their extent is greater than the extent of 
the highest order portion of the system. If their extent is greater, they may require addi- 
tional boundary conditions and this could adversely affect the regularity at the boundary. 
In most problems of interest lower order terms would have an extent no larger than that 
of the highest order terms and the regularity estimates would hold true in the same form 
as (3.3). This is proved in exactly the same fashion as for systems of differential equations. 
That is, the lower order terms can be considered as part of a right-hand side of (2.1) 
and then the estimate (3.3) with this modified data follows. Then using the interpolation 
estimates in Lemma 6.1 the estimate with the original data follows. 

8. Summary of Results and Examples 

In this section we summarize the results of this paper and apply the theory to several 
examples. These examples are chosen to illustrate the theory; the boundary conditions 
we consider are not representative of those used in practice. Unfortunately, more realistic 
boundary conditions lead to a great deal of algebraic manipulation. To apply Theorem 
3.1 to determine the regularity of a boundary value problem for a regular elliptic system 
of difference equations on a boundary-fitted coordinate system one must only consider 
the “frozen coefficient problem” for the system at each boundary grid point. The frozen 
coefficient problem at a point on the boundary is the constant coefficient problem obtained 
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by fixing the coefficients of both the system and the boundary conditions at their values 
at that boundary point. This frozen coefficient problem is considered on the half-space 
determined by the tangent space to the boundary and the inward unit normal at the point. 
If the frozen coefficient problem is regular for each point on the boundary then the original 
variable coefficient problem is regular. Thus we need only consider constant coefficient 
boundary value problems on a half-space. 

The steps one must take to check the regularity of the boundary value problem are 
as follows. The n-tuples a and r defining the order of the elliptic system of difference 
equations must be determined and lower order terms can then be neglected. The regular- 
ity of the scheme must also be checked (Bube and Strikwerda [3]). The reduced equation, 
obtained by Fourier transforming in the tangential variables, must satisfy the resolvent 
condition of Section 2 and must be adjusted so that Assumption 2.1 is satisfied. With the 
reduced equation in this form the number of boundary conditions can be determined (As- 
sumption 2.1) and they should be ordered so that Assumption 2.2 is satisfied. Assumption 
2.3 will be satisfied for most systems arising in practical application. The final step is to 
check for eigensolutions of types I, II and III. Theorem 3.1 states that regularity up to the 
boundary is equivalent to the nonexistence of eigensolutions. 

The regularity estimate (3.3) shows that for those schemes which require as many 
boundary conditions as does the associated differential equation the solution and its finite 
differences are bounded independently of the grid spacing. For those schemes which require 
more boundary conditions, i.e. numerical boundary conditions, the estimate (3.3) shows 
that these should be of high order to achieve smooth solutions. These observations seem 
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to justify the use of compact difference schemes wherever possible. (A compact scheme is 
one which has the smallest extent possible, for a given accuracy.) For those schemes that 
do require numerical boundary conditions, Theorem 3.2 shows that these extra conditions 
should be of sufficiently high order so as not to affect the accuracy of the solution and 
its finite differences near the boundary. The interior regularity estimates of Bramble and 
Hubbard [2] and others show that for second-order elliptic equations the finite differences 
of the solution are approximations of the corresponding derivatives with the same order of 
accuracy as the solution itself away from the boundary. Theorem 3.2 shows that this can 
also be true up to the boundary under certain circumstances. In particular, it can be true 
if no numerical boundary conditions are required and the boundary conditions are of the 
same order of accuracy as the scheme. 

To illustrate the theory consider several examples. We begin with the Cauchy- 
Riemann equations 

u x -v y = fi 

Uy + V x = / 2 , 

on the half-space {(x,y) : x >0, y € R} with boundary condition 

u(0,y) = y(y). 

Define a(T) := (^)(T 2 - AT + 3) and a{T ) := -a(r -1 ). We approximate the elliptic 
system with the second order accurate scheme given by 

(spy spy) («)(*•'-»<•) = f 8 - 1 ) 
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One can easily check that this is a regular elliptic system of difference equations. Let 
f(u/) := a(e tw ) be the symbol of a. Fourier transforming with respect to the tangential 
variable yields the reduced system 


( «(T*) 

WM 



= F{x v ,u), 


v > 2, ui := h £. 

To determine how many boundary conditions are needed by the system (8.1), we must 
consider the eigenvalues of the resolvent equation, (see Assumption 2.1), given by 


det 


( «(*) 

WM 



( 8 . 2 ) 


If u> 7^ 0, Lemma 2.1 implies that none of the four roots of this equation is on the unit 
circle. Since a(z) = — a[z~ l ) we conclude that there are exactly two roots inside the unit 
circle and two outside. Denote the two which are inside by z\ and z%. Because there are 
two roots inside the unit circle, it is necessary to specify two boundary conditions in order 
to have a well— posed system of difference equations. The conditions which we impose are 


a ) |(«o + ui) + (1 - e)u 0 = <7 0> 

b) h~ r {T-\) r v 0 = 0, 


( 8 . 3 ) 


for e real and r a positive integer. The first boundary condition is the operator Si corre- 
sponding to specifying u and the second is S 2 ; we have pi = -1, p 2 = r - 1. 
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Now we consider the resolvent condition. One could take a] 1 " = aj = P 2 = 0,/?* = 
2 , 0 .} = a .2 = P\ = 0,02 = — 2, or one could take af = 2,aJ = = 0, = 

Pi — P 2 = 0 5 a 2 = ~ 2. Following the procedure for reducing to a canonical form, i.e., all 
differences forward, we obtain two equivalent reduced systems which are to be analyzed 
for eigensolutions. We will work with the one which comes from the second set of a and 
p. This system is 


( ~T(0)T 2 a(T)T 2 ) ( 5 ) ^ ~ ( T 2 f 2 ) ^ 

v > 0, 6 := hu. 

If we would have chosen the first possibility then the two variables uo and u i would have 
been declared superfluous. We could then eliminate these by applying the first equation 
at u = 0 and u = 1. 

The general decaying solution to the resolvent equation is 

( “ ) t 1 - “ ( «(*,) ) + C2 ( «(*,) ) < 8 - 4 > 

First we check for eigensolutions of type III. Since one of the roots goes to 1 as u> — *• 0 
the general decaying solution at u = 0 becomes 

*>o. 

This shows that using any extrapolation with r > 1 for a numerical boundary condition 
gives no eigensolutions of type III. Note that if we would have taken (3T — l) ri (T — l) r2 t> 0 = 
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0 for a numerical boundary condition instead of (8.3b) then there would be eigensolutions 
of type III. For example, take rj = 1 and r 2 = 2, then this condition becomes u 0 = 
36 3 — 7v 2 + 5ui and this boundary condition would not be regular. 

Since (8.1) with the indicated boundary condition is regular, there are no eigenso- 
lutions of type II, so it only remains to check for the existence of type I eigensolutions. 
Applying the boundary conditions to (8.4), we obtain the condition for the existence of a 
nontrivial solution as 


det ( £ ( Zl + J )/ 2 + " e ) £ ( Z2 + X )/ 2 + ( J “ £ ) ) = n 

V (*1 - l) r a(zi) (*2 - l) r a(z 2 ) ) 


(8.5) 


Using the two equations a(zj)a(zj) - |f| 2 = 0, for j = 1,2, we obtain the equation 


a(zi)a(zi ) - a(z 2 )a(z 2 ) = 0. (8.6) 

We now have to determine whether there are any roots zi,z 2 of equations (8.5) and (8.6) 
which are both of modulus less than unity. Using the symbolic manipulation language 
MACSYMA [8] to determine the solutions to this equation, we observe that there are no 
solution with z\ and z 2 both less than or equal to 1 in magnitude. Thus we conclude that 
there are no eigensolutions of type I for the cases e 6 {0, l},r € {1,2,3}. Therefore the 
system is regular up to the boundary for these values of e and r. 
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As the next example, consider the Stokes equations 

Au + p x = fi 
Av + p y = f 2 

U X + Vy = / 3 , 

on the half-space {(x, y) : x > 0, y € R} with boundary conditions u(0,y) = <7i(y) and 
v(0, y) = g 2 (y). The scheme we consider is 


/ (^x-f^x- + 6y4_6y_) 

0 

S x -\ 

M 

0 

(6 X + S X - + 5y + 6y_) 

6y- 

V (x„,y M ) = F(x l/ ,y fl ), 

V 6*+ 

by+ 

o J 

[pJ 


u ^ 1, /z 6 Z . 

One can easily check that this is a regular elliptic system of difference equations. Let 
b(T) := (T + T -1 — 2)//i 2 and := ( e — 1 )/h. After Fourier transforming with respect 
to the tangential variable, we obtain the reduced system 


jb(T x ) - |f p O (1 

o Mr.) - kl 2 
V (r« - 1) A s 



I Sl (x„,w) = F(i„,a>), 


( 8 . 7 ) 


u > 1. 


The determinant of the resolvent equation is 


(b{z)~ k| 2 O 

det O b{z) - |f| 2 
\{z-l)/h $ 


(1 - z^)/h\ 
? ) 


= 0, 
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and the eigenvalues are the roots of this equation, 


*± = (2 + A 2 |d 2 ± \/(2 + | 2 ) 2 - 4)/2, (8.8) 

where each root is double. We conclude that there must be two boundary conditions since 
z_ is a double root with modulus less than unity (Assumption 2.1). Consider the following 
boundary conditions: 


£ 

a) ~(uo + u 1 ) + (l-£)u 0 = ffi, 

(8.9) 

£ 

b) 2 ^° *’ 1 ) + (1 — Wo = 92- 

Note that, as u — > 0, z± — > 1 so we have no type III eigensolutions. The general decaying 
solution to the homogeneous difference equation (8.7) is 


v J (xyju;) = ci ( [z- - l)/h J z^ + c 2 ( {{u + 1 )z- - u)/h J z v _ 1 . 


( 8 . 10 ) 


Applying the boundary conditions (8.9), the condition for a nontrivial solution is 


det ( e ^ Z ~ + 1 )/ 2 + f ^ _ o 

det 1 (6(z_ + l)/2 + (1 - 6))(z_ - 1) 6z _ + (1 - 6) ) ~ U - 


After simplification this equation reduces to 


edz 2 _ + 26(2 — e) 2 _ + 4(1 - 6) + e6 = 0, 
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and since z^_ = (2 + /i 2 |f| 2 )z_ — 1 we obtain 


2(1 - s)!s 

2 — e + e(2 + h 2 \$\ 2 ) ‘ 

Noting that z-, as given by (8.8), is strictly positive, we conclude that this equation is not 
satisfied if 0 < e < 1 and 0 < 6 < 1. When 6 = 0 it is easy to see that the determinant 
never vanishes. Therefore there are no eigensolutions of type I. This implies regularity up 
to the boundary for this system of difference equations. If e = 0 and 6 is taken larger than 
(2(\/2 — l)) -1 then there exist type I eigensolutions, so the difference equations are not 
regular up to the boundary in this case. 

Finally, consider the same difference equations but with the boundary conditions for 
the differential system given by 


a) u x (0,y) = 0, 

b) v{0,y)=g{y). 

Approximate these equations by 


a) e(u(xi, y M ) - u(x 0 , y M )) + (1 - e)(u{x 2 ,y fl ) - u(x 0 , y M )) = 0, 
£ 

b ) 2 ( v ( x °’ y ^ + v ( Xl ’ + _ 5 ) w ( x o, y M ) = g(xo, y M )- 

After applying these boundary conditions to the solution (8.10) we obtain 


( 8 . 11 ) 


( £ ( z - ~ !) + (1 “ £ )( z - - 1) e + 2(!-e)z-\ _ n 
V (S(z- + l)/2 + (1 - 6))(z_ - 1) 6z_ + (1 - 6) 


( 8 . 12 ) 
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Taking e = 0 and 5 = 0 gives no eigensolutions of type I. However, if e = 1 and 5 = 0 
then the determinant vanishes for any value of u> 6 R. In this case the constants from 
equation (8.10) can be taken to be ci = 1 and C 2 = 1 — z~. Thus, since there exist type I 
eigensolutions, the scheme is not regular up to the boundary. 
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